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20.  Abstract 

The  first  six  quantities  are  material  variables  which  depend  on  the  type  and 
structure  of  the  elastomers  used.  The  next  three  are  operation  variables  and  the 
last,  geometric. 

H e are  currently  investigating  experimentally  each  of  the  material  variables 
for  a given  number  of  "model"  elastomers  and  are  conducting  heat  generation  ex- 
periments with  these  elastomers  during  cyclic  shear  deformations.  The  experi- 
mental results  will  be  used  to  improve  our  simulation  model  for  heat  generation. 

A second  problem  which  we  feel  is  important  in  connection  with  track  pad  life 
and  performance  is  the  structural  uniformity  of  the  pads  achieved  during  manu- 
facturing. We  are  embarking  on  an  investigation  of  the  problem  of  the  effects  of 
material,  process  and  geometric  variables  during  injection  mold  curing.  Our 
initial  efforts  are  numerical  having  as  a goal  the  prediction  of  crosslink  density 
as  a function  of  space  and  time  during  the  molding  process. 
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ABSTRACT 

One  of  the  major  problems  in  track  pad  performance  and  expected  life 
is  heat  generation  brought  about  by  the  periodic  loading  of  the  pads  and 
their  viscoelastic  (partly  dissipative)  nature.  This  problem  has  been 
simulated  numerically  making  use  of  a number  of  simplifying  assumptions. 
The  results  are  included  In  this  final  report.  Conceptually,  the 
temperature  Increase  can  be  expressed  as 

T a f(  j>  , Cp  , K,  E,  jk,  &E,  o)  , h,  b) 

The  first  six  quantities  are  material  variables  which  depend  on  the  type 
and  structure  of  the  elastomers  used.  The  next  three  are  operation 
variables  and  the  last,  geometric. 

We  are  currently  Investigating  experimentally  each  of  the  material 
variables  for  a given  number  of  "model1'  elastomers  and  are  conducting 
heat  generation  experiments  with  these  elastomers  during  cyclic  shear 
deformations.  The  experimental  results  will  be  used  to  improve  our 
simulation  model  for  heat  generation. 

A second  problem  which  we  feel  is  important  in  connection  with  track 
pad  life  and  performance  is  the  structural  uniformity  of  the  pads  achieved 
during  manufacturing.  We  are  embark  ig  on  an  investigation  of  the  problem 
of  the  effects  of  material,  process  and  geometric  variables  during 
injection  mold  curing.  Our  Initial  efforts  are  numerical  having  as  a goal 
the  prediction  of  crosslink  density  as  a function  of  space  and  time  during 
the  molding  process. 
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A.  HEAT  GENERATION 
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I.  Statement  of  the  actual  problem 


When  a tank  is  in  motion,  the  rubber  track  pads  undergo  shear  and  com- 
pressive deformations  as  they  contact  the  ground  surface  (which  is  either 
smooth  or  rough)  and  as  they  pass  the  region  directly  beneath  the  six 
wheels  of  the  vehicle.  The  compressive  stresses  applied  to  the  track  pads 
could  be  schematically  represented  as  a function  of  time  as  follows: 


TfevioU 
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During  the  ’'deformation  period"  a certain  amount  of  the  mechanical  work 
is  dissipated  into  heat  resulting  in  the  heat  of  the  rubber  pads.  This  is 
because  the  rubber  pads  are  viscoelastic  substances  which  during  a cyclic 
deformation  would  exhibit  hysteresis  shown  below 
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Assuming  that  the  rubbery  track  pads  are  "ideal  rubbers",  then  the  in- 
ternal energy  is  not  a function  of  the  state  of  strain  and  the  first  Law 
after  the  end  of  each  cycle  becomes 

AQ  = AW 


since  dE  = 0 . 

The  thermal  energy  generated  increases  the  average  temperature  of  the 
track  pads.  The  temperature  distribution  depends  on  the  heat  transfer 
conditions  around  the  pad  and  on  the  uniformity  of  the  generated  heat  rate 
by  viscous  dissipation,  (AQ)  above.  A temperature  field  T(xi>t)  will 
result  In  variations  of  the  rubber  pad  properties  with  space  and  time,  and 
may  well  cause  mechanical  deterioration  and  aging  because  of 

(a)  Oxidative  degradation  which  is  accelerated  at  high  temperatures 

(b)  Thermal  degradation 

(c)  Promotion  of  unfinished  cross-linking  reactions,  which  when  they 

occur,  may  render  the  pads  brittle  and  susceptible  to  failure 

(d)  The  increased  probability  of  reactions  with  soil  ingredients 

(e)  The  temperature  dependence,  in  the  absence  of  any  reaction,  of 

the  mechanical  properties  of  the  rubber  pads. 

During  the  "rest  period"  the  rubber  pads  have  an  opportunity  to  trans- 
fer the  heat  generated  in  the  previous  period.  Heat  transfer  is  by  con- 
duction to  the  initial  belt  of  the  track  and  by  forced  convection  past 
their  top  surface. 

In  the  following  we  report  on  our  simulation  work  of  the  above  prob- 
lem. We  simplified  the  physical  situation  so  that  we  end  up  with  a 
tractable  mathematical  problem  the  solution  of  which  gives  useful  guides 
on  the  effects  of  material  and  operation  variables  on  the  heat  generation 
and  heat  transfer  in  rubber  track  pads. 

1 1 . Problem  simulation 

The  following  assumptions  were  made  about  the  physical  conditions  and 
material  response  during  the  actual  track  pad  deformation  history: 

a.  A geometrically  simple  flat  rectangular  body  undergoing  cyclic 
tensile  strains  of  amplitude  and  frequency  is  considered 
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b.  The  mechanical  response  of  the  track  pads  is  assumed  to  be  ade- 
quately described  by  a Voigt-Kelvin  linear  viscoelastic  mechani- 
cal model.  Such  a model  would  be  appropriate  for  small  deforma- 
tions 


Mataved  JaaVOiuetavfc  : 


c.  Only  the  viscosity  of  the  material  (and  that  used  in  the  Voigt 
model)  is  assumed  to  be  temperature  dependent.  All  the  other 
relevant  properties,  p , Cp  , k and  the  heat  transfer  co- 
efficient h are  taken  to  be  temperature  independent.  In 
I reality  they  are  linear  functions  of  temperature,  a mild  func- 

tional dependence  compared  to  the  exponential  temperature  de- 
pendence of  viscosity. 
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d.  The  initial  and  boundary  conditions  for  the  simulation  of  the 
heat  transfer  problems  are: 


i.  oT/az  =0  at  z = 0 - symmetric  temperature 

profi le 


U , 


3T 

az 


+ i 


h(T 


~ V 

■-i  SL 


forced  convective  heat 
transfer  at  both  sur- 
faces 


iii.  T = T at  t < 0 
a 


the  rubber  pads  start 
the  process  at  ambient 
temperature. 


III.  Mathematical  description  of  the  simulation  model 

The  starting  pq^i nt  for  -the  solution  of  any  heat  transfer  problem  is 
the  equation  of  energy,  which  in  terms  of  the  temperature  field  is 

Cv £ , _j>  Cf  v.^t  * k - t( - 5.  ■•'?y 

The  term  on  the  left  is  the  rate  of  accumulation  of  thermal  energy  and 
the  terms  on  the  right-hand  side  are,  respectively,  the  rate  of  heat 
convected  by  flow,  rate  of  heat  conducted,  rate  of  heat  involved  in 
P-V  (expansion)  work,  and  rate  of  heat  dissipated. 

Keeping  the  accumulation,  conduction  and  viscous  dissipation  terms, 
the  only  non-zero  quantities  is  the  simulation  problem,  we  obtain 


01 

Ot 


The  above,  as  well  as  the  boundary  and  initial  conditions  can  be  put  in 
dimensionless  form 
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The  term  -t  * -t  e , responsible  for  the  heat  up,  can  be  cal 

XX  dX  XX 

culated  as  follows:  e = e u cos  wt  and  t = -ye  . Only  the  viscous 

W AA 

contribution  of  the  Voigt-Kelvin  element  is  taken  into  account  because 
the  elastic  part  does  not  contribute  to  the  dissipated  energy.  The  vis 

— KT 

cosity  is  a function  of  temperature  y = a e~  . Letting  yQ  be  the 
viscosity  at  T , we  obtain 
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If  the  viscosity  is  taken  to  be  constant  during  each  cycle,  but  dif- 
ferent from  cycle  to  cycle,  then  the  heat  dissipated  per  cycle  period 


AQ 


2 -beT 
nw  E0  e a 


Choosing  the  dimensionltss  time  interval  At  = --E-  , where 

0 ^ 

At  <<  tp  , assuming  that  the  heal  generated  in  each  period  is  equally 
divided  in  each  dimensionless  time  interval  At  , as  shown 
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W6  obtain  for  each  dimensionless  time  interval  At  , a quantity  of  heat 
generated  (in  dimensionless  form) 


£2  2 -beTa  At_ 


This  problem  can  be  formulated  in  terms  of  finite  differences  (Taylor 
series  expansion)  as  follows: 
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IV,  Results  and  discussion 


The  "design  expression"  for  Heat  generation  is 

2 

AQ  = tt  eQ  u v 

and  one  can  see  the  effects  of  the  material  variable  y and  the  operation 
variables  eQ  , and  u . This  heat  generation  term  is  "in  competition" with 
the  conductive  heat  transfer  rate,  which,  inside  the  pad  is  controlled  by 
the  thermal  conductivity  k and  at  the  pad  surface  by  the  heat  transfer 
coefficient  h . 

The  results  below  indicate  the  effects  of  the  variables  eQ  , u>  , y , 
k and  h on  the  temperature  distribution  after  3,000  cycles.  Specific- 
ally: 

a.  (Fig.  1)  - Keep  yQ  constant,  choose  a small  hQ  value: 

Vary  kQ  to  IQko  and  100kQ  , the  result  shows  very  small 
effects.  That  is, the  heat  generated  can  not  be  dissipated 
to  the  outside  if  ho  is  small. 

b.  (Fig,  2)  - Keep  yQ  constant,  choose  hQ  to  be  nearly 

the  same  magnitude  as  1 00kQ : 

Vary  kQ  to-  10kQ  and  100kQ  . The  temperature  rise  for 

large  kQ  may  decrease  by  about  25%  than  for  small  kQ  . 

c.  (Fig,  3)  - Keep  yQ  constant,  choose  a small  kQ  value: 

Vary  hQ  to  1 0hQ  , and  100hQ  . Because  of  k0  is  small, 

the  heat  still  can  not  be  transferred  outside.  Large  hQ  only 

lowers  the  temperature  near  the  rubber  pad  surface. 

d.  (Fig.  4)  - Keep  yQ  constant,  let  kQ  be  large: 

Vary  hQ  to  10hQ  and  100ho  . The  large  hQ  can  reduce 

the  temperature  rise  in  the  rubber  by  about  25%  than  small  hQ  . 

e.  (Fig.  5,  Fig.  6)  - Keep  kQ  and  hQ  value  constant: 

Vary  viscosity  yQ  to  0.1yQ  ; the  result  shows  a signi- 
ficant effect  on  reducing  the  temperature  rise.  The 


temperature  rise  is  only  about  20%  of  the  original  value,  Of 
course,  larger  kQ  and  ho  will  still  help  in  transferring 
heat  and  lowering  the  temperature  rise. 

f.  Vary  the  frequency  and  the  amplitude  give  the  same  result  as 
varying  viscosity.  It  is  important  to  note  that  AQ  depends  on 
the  square  of  the  deformation  amplitude  eQ  . 

From  the  above  we  can  conclude  that  the  following  are  desirable  design  de- 
cisions on  the  elastomeric  track  pads  and  their  overall  support  system  design: 

a.  Since  increasing  the  thermal  conductivity  reduces  the  temperature 
rise  inside  the  pads,  we  should  investigate  the  effects  of  incor- 
porating conductive  additives  to  the  mechanical  properties  of  elasto- 
mers and  decide  whether  it  is  desirable  to  use  them  or  not. 

b.  Consider  track  fender  (cover)  designs  which  would  increase  the  con- 
vective heat  transfer  during  the  rest  periodt  forced  air  systems 
under  the  fender  may  be  good  design  conditions 


c.  The  role  of  viscosity  is  predominant.  We  should  consider  the  ef- 
fect of  material  parameters, such  as  chain  flexibility  and  cross- 
link density  on  the  value  of  both  the  viscosity  and  the  flow  activa 
tion  energy  b . 

d.  The  effect  of  the  frequency  should  be  related  to  the  field  effect 
of  "rough"  roads,  The  presence  of  obstacles  increase  the  fre- 
quency of  deformation. 

e.  The  large  effect  of  the  amplitude  of  deformation  on  the  rate  of 
heat  generation  may  be  related  to  the  presence  of  sharp  (small 
contact  area)  obstacles.  The  smalior  the  contact  area,  the  larger 
the  pressure  (x  ) and  thus  the  larger  the  deformation.  The 
exact  relationship  between  the  obstacle  sharpness  and  is  not 
clear  because  of  the  viscoelastic  nature  of  the  pad. 


The  above  are  simple,  approximate  design  guides  useful  as  such,  and  in 
providing  us  with  directions  for  extending  the  work  into  the  non-linear  re- 
sponse regime  and  is  devising  the  appropriate  tests  to  compare  with  field 
test  results. 


B.  MOLECULAR  CONSIDERATION  OF  THE  HEAT  GENERATION  PROBLEM 
I . Dependence  of  Elastomer  Properties  on  Molecular  Parameters 

In  addition  to  its  influence  on  chemical  degradation  there  is  an 
additional  effect  of  chemical  structure  on  the  basic  strength  properties 
of  a vulcanized  elastomer,  since  fracture  involves  some  mechanical  rup- 
ture of  molecular  chains  by  rupture  of  chemical  bonds  along  the  chain, 
as  has  been  verified  by  electron  spin  resonance  (ESR)  measurements,  and 
chemical  bond  strength  varies  with  the  type  of  chemical  bond.  One  of  the 
most  important  parameters  of  molecular  structure  is  the  degree  of  cross- 
linking.  The  type  of  cross-linking  is  also  a variable  which  may  be  signif- 
icant in  certain  cases  (e.g.,  where  chemical  degradation  is  involved). 

In  representing  viscoelastic  properties  in  our  matrix,  no  distinction 
is  made  between  linear  and  non-linear  viscoelasticity.  This  is  not  a 
critical  point  where  trends  are  being  discussed,  but  it  immediately  becomes 
critical  v;hen  quantitative  calculations  are  made.  Non-linear  visco- 
elasticity is  one  of  the  difficult  areas  in  the  theory  of  mechanical 
properties  of  polymers.  There  is  no  unified  theory  to  treat  this  type 
of  behavior  (as  there  is  for  linear  viscoelasticity).  There  is  there- 
fore a great  need  for  research  in  this  area,  and  for  the  working  out  of 
even  empirical  methods  to  handle  this  type  of  behavior,  since  it  is  the 
type  of  behavior  encountered  in  many  use  applications  (including  the 
tank  tread  application  under  consideration  here)  and  in  many  of  the  areas 
of  polymer  processing,  such  as  Injection  molding  and  other  fabrication 
methods. 

Another  point  of  importance  is  that  the  matrix  presented  here  refers 
to  the  molecular  parameters  of  the  polymer  itself.  Other  components  are 
usually  also  involved  in  a rubber  compound,  such  as  carbon  black,  antioxi- 
dants, curing  agents,  etc.  However,  since  only  trends  are  specified, 
the  relationships  shown  here  should  remain  valid  even  for  commercial  com- 
pounds containing  other  Ingredients  besides  the  elastomer  Itself. 

In  developing  an  algorithm  to  predict  service  life  from  molecular 
parameters  (an  ambitious,  but  very  attractive  goal,  related  to  the  con- 
cept of  "tailor-making"  molecules),  the  molecular  parameters  that  will 
be  most  useful  are  again  those  referring  to  the  general  characteristics 
of  the  molecules  (the  "physical"  characteristics,  one  might  say)  rather 
than  their  exact  chemical  structure.  Put  another  way,  polymer  chemistry 
Usually  must  be  translated  into  polymer  physics  before  physical  and 
mechanical  properties  can  be  predicted  successfully. 
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In  relating  the  mechanical  behavior  of  elastomers  to  molecular  para- 
meters, six  different  areas  of  description  can  be  distinguished: 

1)  Chemical  structure  (atomic  composition  and  chemical  bonds) 

2)  Molecular  characteristics  (chain  stiffness,  symmetry, 
freedom  of  rotation  of  side  groups,  stereospecificity, 
molecular  weight) 

3)  Intermolecular  forces  (secondary  bonding,  hydrogen  bonding, 
dipole-dipole  Interactions) 

4)  Supermolecular  structure  (cross-linking,  aggregation  and 
crystallization) 

5)  Bulk  properties  (glass  transition  temperature,  thermal 
conductivity,  density,  etc.) 

6)  Mechanical  behavior  (heat  generation,  viscoelastic  proper- 
ties, fracture  strength) 

These  categories  are  sometimes  closely  related,  and  may  even  overlap  in 
some  ways,  but  they  clearly  represent  a hierarchy  connecting  the  basic 
chemical  constitution  of  the  molecule  and  the  mechanical  properties  ob- 
served under  use  conditions.  The  possibility  of  making  this  connection 
is  the  basis  of  the  concept  of  "tailor-making"  polymer  molecules  to  pro- 
duce desired  final  properties.  This  can,  in  fact,  be  carried  out  in 
practice  in  some  cases.  An  example  is  the  successful  development  by 
DuPont  recently  of  Qiana  fiber,  which  was  done  by  a highly  logical 
approach  in  which  the  proper  molecular  parameters  were  brought  together 
in  order  to  produce  a type  of  nylon  fiber  which  was  more  like  silk.  The 
resulting  fiber  obtained  was  exactly  what  was  predicted  and  desired. 

In  an  approximate  way,  the  above  list  of  six  categories  could  be 
reduced  to  two,  in  which  the  first  four  categories  could  all  be  lumped 
as  "molecular  parameters"  and  the  last  two  as  "properties."  Our  Interest 
here  is  basically  in  viscoelastic  (V-E)  properties  and  strength  proper- 
ties, since  these  are  most  closely  related  to  the  problem  of  the  service 
life  of  tank  track  shoes.  However,  heat-up  during  use  depends  not  only 
on  the  viscoelastic  properties  (the  effective  "Internal  viscosity"  used 
In  the  heat  generation  calculations  in  the  other  section  of  this  report) 
but  also  on  such  properties  as  thermal  conductivity  (also  used  in  these 
heat  generation  calculations);  these  other  properties  of  secondary  impor- 
tance must  therefore  also  be  taken  into  account.  However,  as  noted  there, 
the  internal  viscosity  (a  viscoelastic  property)  is  the  property  of  major 
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and  overriding  importance  in  the  heat  generation. 

Present  theory  does  not  allow  a quantitative  calculation  of  proper- 
ties from  molecular  parameters  (this  would  be  an  ultimate  goal  for  poly- 
mer scientists).  However,  solid  and  quantitative  theoretical  connections 
can  be  made  between  the  viscoelastic  properties  of  a material  (as  ex- 
pressed by  the  experimentally  measured  relaxation  time  spectrum)  and  the 
heat  generated  in  a cyclic  deformation  experiment.  The  relaxation  time 
spectrum  Itself,  however,  cannot  at  the  present  time  be  predicted  from 
a knowledge  of  the  chemical  structure  of  the  polymer  molecule. 

Qualitative  connections  can,  however,  be  made  between  molecular 
parameters  (preferably  the  general  characteristics  of  the  molecule 
rather  than  the  exact  chemical  structure)  and  observed  bulk  properties. 

As  a concise  way  of  doing  this,  we  have  set  up  a matrix  relating  signi- 
ficant general  molecular  parameters  to  viscoelastic  and  strength  proper- 
ties, and  also  to  the  other  properties  Involved  In  the  heat  generation 
calculations  given  in  the  other  section  of  this  report.  To  indicate 
trends,  we  have  used  a + sign  to  Indicate  that  the  property  increases 
or  improves  as  the  molecular  parameter  increases,  and  a - sign  to  indicate 
a decrease.  Parentheses  are  used  when  the  effect  is  small,  and  a 0 is 
given  where  there  is  no  appreciable  effect.  This  is  shown  on  Table  I. 

One  important  molecular  parameter  which  cannot  be  Included  in  this 
table  is  the  chemical  composition  of  the  molecule,  since  this  Is  some- 
thing that  cannot  be  described  in  terms  of  plus  or  minus  trends.  However, 
it  Is  a very  Important  parameter  in  terms  of  thermal  and  oxidative  (both 
oxygen  and  ozone)  degradation,  since  these  are  chemical  processes  specif- 
ically related  to  the  chemical  composition  of  the  molecule.  This  deg- 
radation has  a very  serious  effect  on  the  fracture  strength  properties 
of  elastomers;  even  a relatively  small  amount  of  degradation  will  pro- 
duce an  appreciable  decrease  of  tensile  strength.  Here  is  a case  also 
where  the  properties  themselves  Interact,  since  heat  generation  from 
viscoelastic  deformation  will  raise  the  temperature  of  the  elastomer,  and 
this  in  turn  will  speed  up  the  rate  of  chemical  degradation.  Another 
case  of  properties  which  are  related  are  the  glass  transition  temperature 

(T  ) ana  he  V-E  activation  energy;  a lower  value  of  T for  the  polymer 
y y 

will  produce  a lower  value  of  Eacj.  (defined  by  the  Arrhenius  rate  equa- 
tion) for  any  given  experimental  temperature. 


rysta 1 lizatio 
Aggregation 
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II.  Heat  Generation  - Molecular  Phenomena 


In  relating  heat  generation  to  molecular  structure  (which  would  be  one 
of  our  first  goals  In  making  a connection  with  molecular  parameters),  a 
number  of  factors  would  have  to  be  analyzed  and  taken  Into  account.  The 
nature  of  the  chemical  structure  Itself  may  not  allow  a direct  prediction 
of  the  Internal  viscosity  (or  heat  generation)  of  the  rubber.  As  an  Illustration 
of  this  point,  It  can  be  pointed  out  that  two  rubbers  which  are  known  to  show 
opposite  extremes  In  heat  generation  — natural  rubber,  which  has  very  low  heat 
generation,  and  butyl  rubber,  which  shows  unusually  high  heat  generation  — are 
both  hydrocarbon  rubbers.  Our  explanation  which  has  been  offered  for  the  high 
Internal  friction  of  butyl  rubber  Is  the  "bumpiness"  of  the  molecules,  which  con- 
tain very  many  methyl  side  groups,  In  contrast  to  the  natural  rubber  molecule, 
which  contains  relatively  few  methyl  side  groups.  This  suggests  attempting  a 
correlation  with  the  number  and  size  of  the  side  groups  on  the  polymer  chain. 

Some  further  suggestive  evidence  In  this  direction  Is  the  fact  that  SBR  rubber 
(styrene/butad Ine  copolymer)  Is  also  known  to  have  a relatively  high  heat 
generation,  and  this  polymer  has  large  phcrryl  rings  as  side  groups.  Following 
out  this  reasoning,  the  polymer  with  the  lowest  heat  dissipation  would  be 
predicted  to  be  linear  (non-branched)  polyethylene.  Unfortunately  this  polymer 
(usually  known  as  high  density  polyethylene,  or  HD  PE ) Is  cry-.tal  1 Ine,  and  Is 
therefore  a plastic  rather  than  a rubber.  Even  branched  polyethylene  Is  crystalline 
and  not  an  elastomer.  A certain  nuumber  of  small  side  group  seems  to  be  necessary 
on  a polyethylene  chain  to  make  It  a rubber;  copolymerization  of  ethylene  with 
propylene  (In  approximately  50/50  ratio)  gives  a polymer  which  Is  an  elastomer 
(called  EPR  rubber).  A full  analysis  of  heat  generation  from  this  point  of  view 
should  be  carried  out,  and  we  would  plan  to  do  this  In  the  next  phase  of  this 
project. 

The  "bumpiness"  of  the  polymer  molecules  produces  some  "Interlocking" 
of  adjacent  molecules,  which  would  be  an  analog  at  the  molecular  level  "sur- 
face roughness"  which  Is  known  to  be  an  Important  factor  In  the  sliding  friction 
of  solid  surfaces.  This  type  of  Interlocking  mi ght  be  reflected  In  an  Identifiable 
way  In  the  bulk  density  of  the  rubber.  A possible  correlation  of  heat  generation 
with  bulk  density  should  therefore  also  be  checked. 


III.  HEAT  GENERATION  MEASUREMENTS 


(Linear  Sinusoidal  Deformation) 

Equipment  has  been  designed  and  constructed  to  measure  heat  gener- 
ation In  the  linear  range  during  sinusoidal  deformation  of  a rubber  sanple. 

The  Instrument  Is  being  operated  In  shear,  but  can  easily  be  adapted  to 
operate  In  tension,  flexure  or  compression.  A temperature  box  has  also  been 
built  around  the  sample  so  that  experiments  can  be  carried  out  at  different 
ambient  temperatures.  A sketch  of  the  equipment  Is  shown  In  an  attached 
figure. 

The  sample  Is  deformed  by  a 1/8  HP  motor  fitted  with  an  eccentric 
of  adjustable  radius  to  allow  the  amplitude  of  deformation  to  be  varied. 

The  frequency  can  be  varied  by  a variable  speed  drive  attached  to  the  motor 
which  changes  the  RPM  of  the  motor  by  a simple  manual  dial.  The  frequency 
Is  the  same  as  the  motor  RPM,  which  Is  measured  wtth  ± 1%  error  by  an 
accurate  stroboscope. 

The  rubber  sample  used  for  sheer  measurements  Is  In  the  form  of  a 
sandwich  with  two  thick  rubber  blocks  on  opposite  sides  of  a central  driving 
plate  (attached  by  a connect  ng  rod  to  the  motor  eccentric),  with  two 
metal  clamping  plates  on  the  outside  of  the  sandwich.  The  rubber  blocks  are 
one-inch  thick,  made  by  cementing  two  pieces  of  rubber  cut  from  a cured  sheet 
1/2  Inch  thick  with  epoxy  glue.  The  cross-sectional  area  Is  a little  over 
one  spuare  Inch.  The  two  rubber  blocks  are  cemented  to  the  center  plate 
and  the  outside  plates  by  epoxy,  and  a copper-constantan  themocouple  Is 
Inserted  In  the  Interior  of  one  of  the  blocks.  Interval  sample  temperature 
Is  read  using  a leeds  and  Northrup  potentlonmcter . 

Results  of  a heat  generation  measurement  arc  shown  In  an  attached  graph. 
The  exact  quantitative  nature  of  this  graph  should  not  be  overemphasized 
because  this  experiment  was  carried  out  using  a friction  grip  between  the 
rubber  blocks  and  the  central  plate  (with  emery  cloth  Slued  to  the  surfaces 
of  the  central  plate),  and  It  Is  believed  that  some  slippage  took  place  at  thl 
point  which  would  add  a frictional  contribution  to  the  heat  generation,  In 
addition  to  the  Internal  viscoelastic  lossed  (sometimes  called  "internal 
friction"  since  It  arises  from  the  polymer  molecules  rubbing  past  each  other). 
However,  the  general  shape  of  such  a heat  generation  curve  Is  evident.  There 
Is  an  Initial  rise  (almost  linear)  as  heat  starts  to  be  developed,  and  heat 
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loss  from  the  sample  Is  negligible.  There  Is  then  a curving  off  as  heat 
loss  (due  to  thermal  conductivity)  starts  to  become  significant;  and  then  a 
final  leveling  off  of  the  curve  as  a steady  state  Is  reached  in  which  heat 
generation  is  just  balanced  by  heat  loss  from  the  sample.  Under  the  experi- 
mental conditions  used  here,  the  temperature  rise  observed  is  only  a few 
degrees  centigrade.  The  amount  of  heat  generated  could  of  course  be  Increased 
by  using  larger  amplitudes  of  deformation.  This  would  require  the  use  of  a 
heavler-duty  motor  (greater  HP),  since  the  work  of  deformation  would  be 
greater.  This  would  also  take  us  into  the  non-linear  region,  as  we  move 
into  the  region  of  greater  power  and  correspondingly  greater  lossed  (of  work 
or  energy  converted  into  heat).  And  this  is  something  that  we  will,  In  fact, 
want  to  do  as  our  studies  extend  Into  the  non-linear  region.  But  these  linear 
experiments,  at  relatively  small  deformations,  will  allow  comparisons  to  be 
made  with  a linear  viscoelastic  model,  and  related  heat  capacity  and  thennal 
conduct tvety  data. 

The  linear  viscoelastic  mode]  discussed  In  a previous  section  (Voigt 
Model)  Involved  a single  elastic  spring  and  viscous  dashpot  and  can  be 
expressed  in  terms  of  an  effective  viscosity  coeffecient, v , as  a represen- 
tation of  the  "Internal  viscosity"  (molecular  friction  effects)  of  the  system. 

Iri  the  linear  region,  this  value  is  a constant.  Actually,  it  is  known  that  this 
Internal  friction  arises  from  a wholu  spectrum  of  molecular  responses  having 
different  relaxation  times,  A*  The  energy  losses  can  be  expressed  In  a more 
detailed  way  by  a "loss  modulus",  E"„  which  can  In  turn  be  expressed  as  an 
Integral : 

«*  w X 

E"  . P H (\1 dA 

o I+w3  a 

where  the  function  H(X)  is  the  linear  viscoelastic  "distribution  of 
relaxation  times"  function,  w is  circular  frequency  and  \ is  the  relax- 
ation time  (\  In  a simple  model). 

Tills  Is  the  theoretically  complete  and  rigorous  formulation  for  the 
linear  case,  but  In  order  to  use  It,  It  is  obviously  necessary  to  know  the 
distribution  function  H ( A)  • Tills  function  cannot  be  calculated  at  the  present 
time  from  a prior  theory,  but  can  only  be  determined  by  experimental 
measurement.  A necessary  part  of  our  experimental  program  will  therefore 
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be  the  measurement  of  this  relaxation  time  spectrun  for  the  rubber  samples 
used  In  our  heat  generation  measurements.  This  can  be  done  to  some  extent 
using  the  Rhcovlbron  dynamic  mechanical  instrument  which  we  already  have 
(in  which  different  fixed  frequencies  can  be  used,  but  where  a continuous 
temperature  scan  Is  the  primary  variable).  But  this  should  be  supplemented 
by  measurements  In  which  frequency  can  be  varied  as  a primary  variable,  at 
a fixed  temperature.  Although  the  H(X)  function  should  In  principle  be 
known  from  \ “ 0 to  ■*  <»  for  substitution  In  the  above  Integral  , the  multiplying 
function  In  the  Integral  Is  such  that  a knowledge  of  tne  HQ)  function  for  values 
of  (%-  l/u>)  one  or  two  orders  of  magnitude  above  and  below  the  experimental 
frequency,  , Is  all  that  Is  necessary  in  actual  practice. 

We  will  eventually  wont  to  measure  hon-1 Ineor  viscoelastic  spectra 
as  a function  of  frequency  and  temperature.  Commercial  Instruments  ore  now 
available  which  are  able  to  do  such  measurements. 
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C . Mathematical  Fitting  of  Non-Linear  Stress-Strain  Curves 

♦ 

One  approach  which  is  being  developed  to  treat  the  non-linear  case 
is  an  empirical  mathematical  fitting  of  the  stress-strain  curve  using  a 
power  series.  By  fitting  both  the  extension  and  retraction  parts  of  the 
stress-strain  cycle  by  such  a series,  the  energy  loss  (or  heat  generation) 
per  cycle  can  be  calculated  as  the  difference  in  area  under  the  two  curves, 
and  this  can  also  be  expressed  as  a power  series,  where  the  coefficients 
are  simply  obtained  by  subtraction  of  the  corresponding  coefficients  for 
the  extension  and  retraction  parts  of  the  cycle. 

This  is  a very  simple  and  convenient  type  of  mathematical  representa- 
tion, and  the  question  is  how  adequately  actual  experimental  stress-strain 
cycles  can  be  fitted  in  this  way.  We  plan  to  measure  stress-strain  cycles 
ourselves  on  the  Instron  machine,  as  soon  as  rubber  samples  become  avail- 
able. However,  we  have  made  some  preliminary  tests  of  the  method  using 
stress-strain  cycles  from  the  published  literature.  It  is  known  that  the 
stress-strain  relation  for  an  elastomer  is  very  non-linear,  with  a strongly 
pronounced  S-shape  when  the  stress  or  strain  is  taken  to  high  values. 

The  non-linear  shape  at  lower  stresses  and  strains  can  be  interpreted  by 
the  Theory  of  Ideal  Rubber  Elasticity  (developed  by  Ruhn  in  the  1930's). 
This  predicts  a stress-strain  function  of  the  form: 

o = N k T [a  — 
a 

where  o is  stress,  N is  the  number  of  "network  chains"  per  unit  volume 
(these  are  the  sections  of  molecules  between  cross-links),  k is  Boltzmann's 
constant,  T is  absolute  temperature  and  a is  the  extension  ratio,  or  ratio 
between  stretched  and  unstretched  length  of  the  sample  (a  = ^ A0 ) • This 
theory  applies  to  a vulcanized  rubber,  and  describes  the  equilibrium  elas- 
tic behavior  only,  not  the  viscoelastic  behavior.  This  therefore  gives 
no  description  of  the  heat  generation  phenomenon,  since  the  extension  and 
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and  retraction  parts  of  a stress-strain  cycle  would  be  identical,  accord- 

» 

ing  to  this  theory,  so  that  there  would  be  no  energy  loss  or  heat  genera- 
tion. 

A stress-strain  function  of  slightly  different  shape  is  obtained  by 
5 

modifying  the  above  equation  into  the  form: 

a C„ 

r‘c,  +~ 

[a  " —J-]  a 

a 

This  is  known  as  the  Mooney-Rivl  in  Equation.  However  this  still  repre- 
sents only  the  first  curvature  of  the  overall  S-shape.  The  large  up- 
swing of  the  curve  at  high  strains  is  usually  attributed  to  a combination 
of  two  effects:  crystallization  of  the  rubber,  and  the  fact  that  the 
network  chains  are  reaching  their  limit  of  extensibility.  The  relative 
importance  of  these  two  factors  will  be  different  for  different  rubbers, 
since  the  tendency  to  crystallize  varies  considerably  for  different  rubbers. 
No  quantitative  theory  has  been  proposed  to  explain  this  final  up-moving 
region  of  the  stress-strain  curve. 

It  is  clear,  in  addition,  that  no  elastic  theory  (even  if  one  were 
available  which  could  represent  the  entire  S-shaped  stress-strain  curve) 
would  be  able  to  predict  the  heat  generation  effect.  For  this,  a differ- 
ence in  extension  and  retraction  curves  is  necessary,  and  this  type  of 
behavior  requires  a viscoelastic  theory.  A viscoelastic  theory  for  this 
high-strain  behavior  is  similarly  not  available  at  the  present  time.  In 
view  of  these  various  difficulties  in  applying  currently  available  rubber 
elasticity  theory  to  the  problem  of  heat  generation  during  cycling  to 
large  strains,  the  idea  of  developing  a purely  empirical  approach  which 
provides  the  needed  mathematical  representation  and  which  can  then  be 
used  in  practical  calculations  seems  to  be  a very  reasonable  and  useful 
t pproach. 


In  checking  the  effectiveness  of  the  power  series  representation  of 
actual  high-strain  experimental  curves  (where  the  S-shape  is  very  pro- 
nounced), a power  series  representation  up  to  the  sixth  power: 

2 3 4 5 6 

a = aQ  + a-j  x + a2  x + a3  x + x + ag  x + ag  x 

was  first  tried,  and  the  accuracy  of  representation  was  found  to  be  sur- 
prisingly good.  The  term  aQ  is  always  equal  to  zero  and  can  be  dropped, 
since  the  curve  starts  at  the  origin.  The  coefficients  a.,  are  obtained 
by  choosing  six  points  on  the  experimental  curve  to  be  fitted  exactly, 
and  then  solving  the  corresponding  matrix  with  a computer  program.  This 
calculation  can  be  carried  out  so  easily  that  the  polynomial  can  be  readily 
extended  to  additional  terms.  We  have  also  checked  the  fit  obtained  using 

g 

a polynomial  Up  to  the  eighth  power  (a Q x ) and  up  to  the  tenth  power 
(a^Q  x10).  The  fit  is  a little  improved  by  using  the  eighth-power  rather 
than  the  sixth-power  polynomial,  but  the  use  of  a tenth-power  polynomial 
seems  to  give  very  little  additional  improvement.  The  eighth  power  poly- 
nomial gives  a very  satisfactory  fit. 

The  sort  of  fit  obtained  is  illustrated  in  the  attached  figure  show- 
ing the  first  stress-strain  cycle  of  a natural  rubber  sample  (taken  from 
Harwood,  Mullins,  and  Payne,  Rubber  Chem.  Tech.,  39,,  814  (1966)),  using 
an  eighth  order  polynomial.  The  strain  parameter  used  was  the  fraction 
of  the  maximum  strain  (x  = e/e  ) in  this  particular  case.  The  points 
used  for  curve  fitting  are  shown  as  circles,  and  the  solid  curves  are  both 
experimental  and  calculated  curves,  since  there  is  essentially  a negli- 
gible difference  between  them  (experimental  and  calculated  curves  super- 
impose). The  coefficients  obtained  for  the  extension  and  retraction  parts 
of  the  cycle  are  listed  below  the  graph. 
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A similar  fit  for  the  second  cycle  of  extension  (not  shown  here)  was 
equally  satisfactory.  There  is  a gradual  change  of  shape  of  the  cycle  as 
cycling  continues,  but  this  finally  settles  down  to  a steady-state  shape, 
corresponding  to  the  experimental  conditions  (here,  a constant  rate  of 
length  change,  as  in  most  stress-strain  tests  today).  This  steady  state 
corresponds  to  a constant  loss  per  cycle. 

Further  investigation  of  this  method  of  curve-fitting  is  being  carried 
out.  The  shape  of  the  stress-strain  curve  is  affected  by  strain  rate, 
temperature,  degree  of  maximum  strain,  elastomer  type,  and  degree  of  cross- 
linking.  The  usefulness  of  this  power  series  method  can  be  extended  by 
learning  how  to  take  these  factors  into  account.  The  results  obtained 
thus  far  verify  that  this  approarh  can  be  a very  useful  one,  particularly 
in  handling  realistic,  high-strain  situations  such  as  are  encountered 
with  tank  treads  and  with  other  heavy-duty  service  conditions. 
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D.  EFFECTS  OF  MANUFACTURING  ON  PAD  STRUCTURE 

During  the  injection  molding  of  track  pads,  a prepolymer,  or 
linear  polymer,  together  with  additives  and  a crosslinking  agent  are 
loaded  into  hot  molds  where  "curing"  occurs.  The  reacting  pad  is  under 
non-uniform  temperature  during  molding  because  of  the  heat  transfer  and 
the  exothermic  heat  of  reaction.  This  fact  leads  to  the  real  possibility 
of  structural  (crosslink  density)  non-uniformities  of  the  pads.  These 
non-uniformities,  in  turn,  may  contribute  to  the  pad  performance  and 
failure. 

The  problem  may  be  simulated  by  considering  a mold  very  wide  in  the 
x and  y directions  and  of  thickness  2H  in  the  z direction.  Assuming: 

a.  A well  mixed  heating  fluid 

b.  No  flow 

c.  No  diffusion 

d.  Constant  conductivity, 

we  arrive  at  the  following  equations  governing  the  reactive  species  and 
energy  balances. 

dc/i 

it'  -ra 

fcP  ^ + *Wa 

n 

where  R^  = Kr  cA  , the  reaction  rate 

= Is  the  heat  of  raction 

K„  = K°e  ^r/RT  is  the  reaction  rate  constant. 
r r 


The  initial  conditions  are 


T(x1 , 0)  * Tg 

CA(xi * 05  = CAq 

The  thermal  boundary  conditions  are: 

^x,  y,  z * 0,  t)  - 0 (symmetry) 
and 

T(x,  y,  t,  z = *H)  = Tw  (isothermal) 
or 

i (x,  y,  t,  z = iH)  = 0 (adiabatic) 
oz 
or 

y.  t,  z . *H)  ■ h(T2=fH- Tnold) 

where  h * is  the  heat  transfer  coefficient. 

Experimentally,  we  need  to  determine  K,  Cp  and  H and  RA-  The 
last  two  quantities  we  plan  to  determine  using  the  Differential  Scanning 
Calorimeter  and/or  the  Mooney  Viscosimeter. 

The  goal  of  this  simultaneous  work  is  to  yield  a relation  such  as: 

CA(x1,  t)  - J (x1?  t)  = f (K,  Cp,  Kr,  n,  Hr,  Tw,  h) 
where:  V*  crosslink  density. 

which  would  be  the  "molding  design  equation."  Work  is  in  progress  in  this 


area. 
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APPEFDIX:  Values  of  material  and  operation  variables  used 
in  the  simulation;  dimensionless  vai'iables 


VIS:  viscosity 
CP:  heat  capacity 
DEN;  density 

CON: thermal  conductivity 


- \e. 


4 /ter  c kk 


C'v  - - ■ s C<x'  /w-^ 

f°  * . v<^ 

*r  S y \u 


CVE:  convective  heat  transfer  coefficient 


modulus 

AE:  activation  energy 

AMP:  amplitude 

FRE:  frequency 

TA:  temp,  of  surrounding 

TING:  increment  time 

TFIN:  final  time 

DIS:  length  of  the  sample/2 

PER:  period  of  oscillation 


k 

£ * o.o  r 

K)  ' >»o 

Tu  - V5T  °C 

\ "Ce  c 1 4 >» 

Jrs  * -»Ows 
-t  ^ * o.o  \ -Se  c.. 


U.  = $ sbw  Vo  U 


Dimensionless  group* j 
Viscosity:  Jj!  - 

Thormal  conductivity: 


/ _ 


VC. 

Convective  coefficient;  L'  - 

U ‘ k. 


Temp. : 

Time: 

Length: 


fe  « 


„ Jtk.  (*  s 


PC 


^ - 
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